iV-particle Bogoliubov vacuum state 
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We consider the Bogoliubov vacuum state in the number-conserving Bogoliubov theory proposed 
by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. We show that in the particle representation 
the vacuum can be written in a simple diagonal form. The vacuum state can describe the stationary 
TV-particle ground state of a condensate in a trap, but it can also represent a dynamical state when, 
for example, a Bose-Einstein condensate initially prepared in the stationary ground state is subject 
to a time-dependent perturbation. In both cases the diagonal form of the Bogoliubov vacuum can 
be obtained by basically diagonalizing the reduced single particle density matrix of the vacuum. We 
compare TV-body states obtained within the Bogoliubov theory with the exact ground states in a 
3-site Bose-Hubbard model. In this example, the Bogoliubov theory fails to accurately describe the 
stationary ground state in the limit when TV — > oo but a small fraction of depleted particles is kept 
constant. 



I. INTRODUCTION 

Bose-Einstein condensate (BEC) is a state of TV bosons 
where all particles occupy the same single particle state 
PJ . Such a state can be an eigenstate of the system only 
in an ideal case when particles do not interact. However, 
dilute atomic gases offer a possibility for experimental 
realization of states that are not too far from ideal BEC 
0- In a dilute system the effect of interactions is weak 
and resulting quantum depletion of a condensate is not 
greater than 1% in a typical experiment [3j. 

From the point of view of theoretical description of 
a system weak interactions allow one to employ various 
perturbation approaches. The approach commonly used 
in the field is the Bogoliubov theory which basically re- 
lies on substitution of a full quantum many-body Hamil- 
tonian by an effective Hamiltonian which is a quadratic 
approximation of the exact Hamiltonian 0, . In the 
original version of the Bogoliubov theory a condensate 
state is identified with a coherent state that for a system 
of massive particles, in non-relativistic quantum mechan- 
ics, can not be acceptable. In the literature one may find 
several versions of number conserving Bogoliubov theory 
[E IE IE IE HI • I n the present publication we consider the 
theory proposed by Castin and Dum [E ■ In this intro- 
ductory part we will briefly describe the key elements of 
the theory. 

Let us begin with a case where N particles, interacting 
via zero range potential V{r) = go6(f), are trapped in a 
time-independent potential U(r). We would like to get 
approximation for the ground state of the system. In the 
Castin and Dum theory bosonic field operator ip is de- 
composed into an operator ao that annihilates a particle 
in a condensate mode cpo and an annihilation operator 
(5-0 defined in the subspace orthogonal to (po, 



ip(r) = (p {r)a + 5ip(r). 



(1) 



We are interested in TV-particle states where a contribu- 
tion coming from the subspace orthogonal to <po is small, 



d 3 r(<5^(r) t 1 5^(r)) = dJV, 



(2) 



where dTV/TV <C 1. Substituting QJ into the Hamilto- 
nian of a system and performing expansion in a "small 
parameter Sip" up to second order leads to the stationary 
Gross-Pitaevskii equation for (po 



H, 



GP<T>0 



2m 



U(r)+g N\fo\ 2 -fi 



0, 



(3) 

and to an effective Hamiltonian quadratic in Sip. Castin 
and Dum have collected operators dtp and ao in such a 
way that the effective Hamiltonian, 



H cS = E (N) + ^J d 3 r (At, - A) £ 
depends only on operators 

a f 

A(r) = -%6iP(r). 



A 

At 



(4) 



(5) 



To get (£fj, sometimes it was necessary to introduce by 
hand the operator ao/y/N, i.e. a quantity of the order 
of one, see (J2J. Diagonalization of Q reduces itself to 
diagonalization of a non-hermitian operator 



C 



H GP + g NQ\cp \ 2 Q 



goNQcftt 



-g NQ*cP* 2 Q -H^ p -g a NQ*\<p \ 2 Q* 



(6) 

where Q = 1 — \<po)(cpo\. Generally the L operator is 
diagonalizable and, for the ground state of the system we 
are interested in, its eigenvalues E m are real. Owing to 
the symmetries of C (i.e. a x Ccr x = — C* and a z Ca z = C\ 
where <j x , <j z are Pauli matrices) it is not difficult to show 
that if 



K) = 



(7) 



is a right eigenvector of C corresponding to an eigenvalue 
E m , the left eigenvector is 



(a a ) 



TV. 



(8) 
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and 



cr x \ip- 



R\* 



(9) 



is a right eigenvector corresponding to —E m . These prop- 
erties imply that the eigenvectors can be divided into two 
classes. The first class (the so-called family "+") consists 
of eigenvectors which can be normalized so that 



(V^M^m) = (u n \u m ) ~ (v„\v m ) = S n 



(10) 



Employing transformation @ we obtain, from eigenvec- 
tors belonging to the family " +" , eigenvectors that form 
the family "-" where (ipf\a z \ipf) -- 
always two eigenvectors of L 



Sij. There are 







corresponding to zero eigenvalue. Collecting all eig> 
tors of C one can decompose the identity operator 
following form 



(11) 

envec- 
in the 



1 = 







«<A) |,0) 







E 

mG"+" 

K 

\u* 



\u m ) 

\v m ) 



{{U m \, 



(o, <^l) 

v m \) 



From the analysis performed it is clear that 

(■u m \<j> ) = 0, 

(VmWo) = 0, 



(12) 



(13) 



and in the subspace of functions orthogonal to a conden- 
sate wave function (j)Q the completeness relation can be 
written in the following form 

1= Yl (\ u ™)(u m \- K)«„|). (14) 

me"+" 

Having calculated Bogoliubov modes (u m , v m ) one can 
easily find diagonal form of the effective Hamiltonian 

Eq(N)+ Emb^bm, (15) 

me"+" 

where so called quasi-particle annihilation operators 

b m = — != (&l(u m \Sip) - a Q (v m \Si/j^)) , (16) 



(17) 



fulfill a bosonic commutation relation 



[b m ,bl 



ml ~ u nm 



■0(1/N). 



The ground state of the Hamiltonian (|15|l is a Bogoliubov 
vacuum state that is annihilated by all quasi-particle an- 
nihilation operators 



b m \0) B =0. 



(18) 



Excited states can be generated by means of the quasi- 
particle creation operators 



|mi,TO 2 , 




(19) 



k=l 



The advantage of the number conserving version of the 
Bogoliubov theory with respect to the original theory is 
that one always deals with states of exactly N particles 
— application of quasi-particle operators does not change 
number of particles thanks to the presence of the So and 
a\ operators in Eq. l|16|) . 

Having a condensate in the ground state of the trap- 
ping potential U (r) (i.e. in the Bogoliubov vacuum 
state |0)b) we may apply a time-dependent perturbation. 
Then 4>q starts evolving according to the time-dependent 
Gross-Pitaevskii equation and the Bogoliubov modes 
evolve according to the time-dependent Bogoliubov-de 
Gennes equations 



ih- 



Uk{r,t) 
dt V v k (r,t) 



C(t) 



Uk(f,t) 

Vk(r,t) 



(20) 



However, the system remains in a Bogoliubov vacuum 
state but this Bogoliubov vacuum and the quasi-particle 
operators b m (t) depend on time. 

Actually one may define any initial Bogoliubov vacuum 
state (not necessary the solution for the ground state of a 
time- independent problem). To this end one has to define 
(f>o and initial Bogoliubov modes (u rn , v rn ) in such a way 
that (|10f> and (|13f) are fulfilled. For example choosing 
4>a = "00 and {u rn — ip m ,v m = 0) for the family "+" 
(where ipj are the eigenstates of the harmonic oscillator) 
the initial Bogoliubov vacuum corresponds to a perfect 
condensate where all particles occupy the ground state 
wavefunction of the harmonic oscillator. 

We should mention that in practice the projection op- 
erator Q in the C operator © can be omitted. One can 
solve the Bogoliubov-de Gennes equations in the stan- 
dard form and at the end apply the operator Q 0, § . 



II. RESULTS 

In this section we will restrict ourselves to analysis of 
the Bogoliubov vacuum state. The Bogoliubov vacuum 
is a state of N particle system where there is no quasi- 
particle, i.e. it is annihilated by all quasi-particle anni- 
hilation operators b m (quasi-particle vacuum), see i|18[) . 
The particle vacuum (a state with no particles) is also 
annihilated by all b m because of the presence of a parti- 
cle annihilation operator in each part of the b m operator, 
see (111)1) . We have shown j^, Hfl DjJ that the Bogoliubov 
vacuum can be obtained from the particle vacuum with 
the help of some particle creation operator d\ 



|0)i 



</ f ) M |o). 



(21) 
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Indeed, if we require that the i$ commutes with all quasi- 
particle operators, 



[b m ,tf]=0, 



then 



b m (#) M |0) = (#) M S m |0) = 0. (23) 
It turns out that two-particle creation operator 



m . 



a, (3=1 



(24) 



(where a^'s create atoms in modes <p a 's orthogonal to 
the condensate wavefunction <f> ) solves the set of equa- 
tions Substituting J23J) to yields the following 
equation for the Z a p matrix 



(3=1 



(25) 



One may wonder how it is possible that a bosonic 
quasi-particle annihilation operator may annihilate, in 
fact, an infinite number of states (for each N there is 
a Bogoliubov vacuum) 12] . It is due to the fact that, 
strictly speaking, b m are not bosonic operators. How- 
ever, for large N corrections to bosonic commutation re- 
lations are small [see Q17fl] and at the present order of 
approximation we may consider b m as bosonic operators 
that allows us to define (within a subspace of a fixed N) 
a Fock space (fB3|) . 

Let us now show that solution of l|25(l can be signifi- 
cantly simplified thanks to properties of the Bogoliubov 
modes 11]. To this end let us multiply l(2*5*|) by (ip*\v m ) 
and then sum over m which leads to the following equa- 
tion 

oo 

( v ;\dp\< P * a ) = Y i { v ;\A\ v >p)Zp a , (26) 

(3=1 

where dp and A operators are defined as 

dp = \v m )(v m \, 
m£"+" 

A = Yl M<«m|. (27) 



mG"+" 



Single particle density matrix corresponding to the Bo- 
goliubov vacuum state reads 

p(f,F')=Nati(?)<l> (f')+ E '), (28) 

77lG"+" 

which indicates that dp operator is a part of the single 
particle density operator corresponding to the subspace 



orthogonal to the condensate wavefunction </>q. Thus 
eigenstates of the single particle density matrix 



(22) p (f } f) = No ^(r)0 o (r") + E diV « KWM?'), (29) 



(where dN a are numbers of particles depleted from a con- 
densate) diagonalize also dp 



dp=Y,dN a \<j>* a ){K 



(30) 



Employing the completeness relation (|14|) and the or- 
thogonality relation l|l(J|) it is easy to show that 



dp A = Adp*, 



(31) 



which in turn implies that eigenstates of the single par- 
ticle density matrix allows us to find also diagonal form 
of the A operator 



a = £a q |O(0 q | 



(32) 



a = l 



There is another equality fulfilled by the dp and A oper- 
ators 



A* A = dp* (1 + dp*), 



(33) 



which allows us to get simple relation between absolute 
value of Aq, and a number of particles dN a depleted from 
a condensate: 

|A Q | 2 = dN a {l + dN a ). (34) 

Let us turn back to the equation (|26() . If we choose as 
the basis vectors (p a the eigenstates 4> a of the single par- 
ticle density matrix we will obtain dp and A in diagonal 
forms and consequently also Z a p matrix in a diagonal 
form, 



Z a (3 — A a <5 Q( 3, 



where 



Note that 



\K\ = 



dN Q 



dN a 



(35) 



(36) 



(37) 



dN a + l 

Finally the Bogoliubov vacuum state in the particle rep- 
resentation can be written in a nice diagonal form 



N/2 



|0)i 



10}, 



(38) 



where aj, operators create particles in modes that are 
eigenstates of the single particle density matrix. 

The analysis performed shows beautiful properties of 
the Bogoliubov modes. Almost all we need in order to 
obtain the Bogoliubov vacuum is given by eigenstates 
and eigenvalues of the single particle density matrix. The 
only things that can not be determined from this matrix 
are the phases of A a . 



4 




N number of a state 



FIG. 1: Panel (a) shows number of atoms depleted from a 
condensate wavefunction, i.e. diV, versus N for Q — 0.017V, 
calculated in the exact state (full circles), calculated directly 
in the Bogoliubov vacuum state in the particle representation 
(open circles) and calculated with the help of Eq. 1451 (open 
squares). In panel (b) the squared overlap between the exact 
ground state of the system \ipo) and the Bogoliubov vacuum 
|0)b in the particle representation is presented. 



III. COMPARISON OF THE BOGOLIUBOV 
THEORY WITH EXACT SOLUTIONS FOR A 
MODEL SYSTEM 

Knowledge of the Bogoliubov vacuum state in the par- 
ticle representation allows for precise comparison of the 
Bogoliubov theory with exact solutions for model sys- 
tems. In Ref.|9J we have analyzed N particles in a dou- 
ble well potential within the tight binding approximation 
where we have found that, quite surprisingly, the Bo- 
goliubov theory gives extremely good predictions for the 
ground state and low lying excited states in the entire 
range of the system parameters. This is rather an excep- 
tion than a rule and therefore in order to test the theory 
we will focus now on a slightly more complex system, 
i.e. N particles in a triple well potential. Restricting the 
single particle Hilbert space to the ground states in each 
well only, the Hamiltonian of the system reads 



H 



(id) 



\A - 1), (39) 



where Cj denotes bosonic operator that annihilates a par- 
ticle in an i-well, VL stands for tunneling rate between 
neighboring wells. This is the only parameter because the 
particle interaction coefficient was chosen as energy scale. 
For VL ^> N the ground state of the system is a BEC 
where all atoms occupy the same condensate wavefunc- 
tion On the other hand Q < 1/N defines the Mott 
insulator regime where the Bogoliubov theory should def- 
initely stop working. 



FIG. 2: Panel (a): the solid line corresponds to the exact 
ground state |^o), the dashed line to the Bogoliubov predic- 
tion |0)b and the dotted line to the perfect condensate state 
I TV : 0o ) where all atoms occupy the condensate mode 140L 
for N = 20 and Q = 1 (dN = 0.4, |(V>o|0)b| 2 = 0.99 and 
|{^>o|iV : 0o)| 2 = 0.83). Panel (b): the solid line corresponds 
to the exact ground state \ipo) and the dashed line to the Bo- 
goliubov prediction |0) B for N = 20 and Q = 0.1 (dN = 2.5, 
|(^o|0)b| 2 = 0.76). The states are projected on |ni, 712,713) 
where m is a number of particles in an i-well. 



The ground state of the Gross-Pitaevskii equation cor- 
responding to the triple-well system reads 



Ao = ^= (1,1,1) 



(40) 



and the solutions of the Bogoliubov-de Gennes equations 
arc 9] 



X(f>± -<f>± 
u± = , v± = 



VX 2 - 1 



where 



Fx V 



V3 



y/X 2 - 1 



/3 e - 4 27r/3 



(41) 



0_ = i=(l, e ^/ 3 ,e^/ 3 ) 



(42) 



and 




Diagonal form of the Bogoliubov vacuum in the particle 
representation is defined by the d) operator 



1 „ t -t 1 
a 6«6 ~ Y a i a i + ^ a 2 a 2: 



(44) 



5 



0.06 




FIG. 3: Panel (a) shows fraction of depleted atoms from a 
condensate wavefunction, i.e. dN/N, versus N for Q, chosen 
so that dN/N w 2/[N(X 2 - 1)] = 0.05, calculated in the 
exact state (full circles) , calculated directly in the Bogoliubov 
vacuum state in the particle representation (open circles) and 
calculated with the help of Eq. 1451 (open squares). In panel 
(b) the squared overlap between the exact ground state of the 
system \tpo) and the Bogoliubov vacuum |0)b in the particle 
representation is presented. 



of f2 are shown in Fig. When is sufficiently large we 
get perfect agreement between the Bogoliubov and exact 
results. 

It is generally believed that when the fraction of atoms 
depleted from a condensate is small, say a few percent, 
the Bogoliubov theory is reliable. Let us now investigate 
the limit when N increases but dN/N remains constant. 
In Fig. [3] we present the results for different N but for 
dN/N w 2/[N(X 2 - 1)] = 0.05. One can see that the 
greater N the worse the Bogoliubov predictions. Actu- 
ally increasing N and keeping dN/N fixed we go far and 
far away from the BEC regime where fl^> N. Indeed, 




O (1/N) , 



(46) 



and in order to keep dN/N constant the f2 must behave 
like 1/N. This example shows that even if the fraction of 
depleted atoms is very small and single particle density 
matrix can be estimated quite satisfactory, for very large 
number of atoms the Bogoliubov theory fails to predict 
structure of the TV-body ground state. Basing on the 
triple-well problem only it is dificult to judge how general 
is such a behaviour. 



where al creates a particle in the condensate mode 



while at and 



o 

2 create particles in 4>i = (<fi + +</>_)/ y/2 
and 4>2 — (<f>+ — 0_)/v2 modes, respectively Num- 
ber of atoms depleted from the condensate mode can be 
directly calculated from the Bogoliubov vacuum in the 
particle representation but it is given approximately by 



dN 



dx (S^Stp) 



dx (Sffiao —OqS'P) 



dx (A 1 A) = (v+\v + ) + (v-\v-) 
2 



X- 



1 



(45) 



Increasing N but keeping fixed number of atoms de- 
pleted from the condensate mode one approaches the 
limit where the Bogoliubov theory should become more 
and more accurate. We have done calculations for differ- 
ent N but for fixed Q/N chosen so that dN w 2/(X 2 - 
1) = 1.5. The squared overlap between the Bogoliubov 
vacuum state and the exact ground state of the Hamil- 
tonian (|39[) is depicted in Fig. ^ where one can indeed 
observe that with increasing N the Bogoliubov solution 
approaches the exact one. With increasing N also num- 
ber of depleted atoms calculated directly from the Bogoli- 
ubov vacuum state and with the help of the estimation 
(I45|l coincide with the exact results. Example of the Bo- 
goliubov and exact states for N — 20 an for two values 



IV. CONCLUSIONS 

We have considered the Bogoliubov vacuum states in 
the number-conserving version of the Bogoliubov theory. 
It is shown that the vacuum can be obtained in a simple 
diagonal form in the particle representation. To this end 
one has to basically diagonalize single particle density 
matrix of a system only. The Bogoliubov vacuum can 
describe eigenstates of a time independent system and 
also a dynamical case where starting with a system in the 
ground state particles are perturbed by time dependent 
force. 

TV-body states obtained with the help of the Bogoli- 
ubov theory have been confronted with exact solutions 
for a triple-well system. In the limit of N — > oo and 
dN = const., the Bogoliubov vacuum states approach 
the exact solutions as expected. However, it turns out 
that in the limit of N — -> oo and for arbitrary small but 
fixed dN/N, the Bogoliubov theory is not able to predict 
iV-body states of the system. 
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